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The aim of this report is to study the behavior of some statistical features of currency rates 

dynamics. We shall also analyze the applicability of the celebrated Samuelson model and its 

modification with normal inverse Gaussian (NIG) increments of logarithms to currency rates 

,— i \ dynamics modeling. We shall always work not with currency rates themselves, but with their 

r"* ' logarithms increments. 

We took rates of 12 currencies from the website [6]. These are the dollar, the euro, the Australian 

,—t, ! dollar, the Canadian dollar, the Iceland krona, the Kazakhstan tenge, the Norway krona, the SDR, 

the Singapore dollar, the UK pound, the Ukranian griven, the Japanese yen. The rates ar taken 

with respect to Russian rouble for the period from January 1, 2000, to December 31, 2009. One 

C\| ■ value of each currency rate is taken for every working day. 

Each of these rates is analyzed for each of the 20 half-years which comprise this 10-year period, 
f— i \ In addition, similar analysis is carried out separately for the period before the financial crisis and 
the period during this crisis. We conventionally consider August 1, 2008 as the beginning of the 
^ . crisis. 

The simplest model for currency rate is the Samuelson model X{t) = e at+aW ^\ where W is the 
standard Brownian motion. But here we operate with discrete time and hence we deal with the 
sequence Xq, . . . , X/y, where Xk = ln(Xk+i/Xk), h — 0, . . . , N — 1 is a selection from the normal 
i-H > distribtion A/"(a; cr 2 ). 

1. Empirical Volatility. 
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This s 2 dramatically rose during the crisis (its order substantially increased). One may consider 
the volatility as an indicator of weight of the financial situation. It is noteworthy that different 
currencies have very different s 2 (with different orders). 

This Samuelson model has a drawback, namely that Xk are Gaussian. Practice shows that the 
tails of these variables have heavier tails. We shall test whether these variables are Gaussian. One 



could estimate mean and variance and apply the standard Kolmogorov test. But we shall perform 
- ■ other tests. 

2. Empirical Asymmetry and Excess Coefficients. The results of this test show a relatively 
poor conformation with the assumption that Xt are normal. These coefficients show the largest 
deviation from the theoretical values for the dollar during the pre-crisis period and for the 
Kazakhstan tenge during the crisis period. 

3. Pickands Estimate. (Cf. [5], chapter 3.) This is an estimate of the extreme value index 
of a selection distribution. This index is equal to for the normal distribution. These results do 
not conform with the Samuelson model too, and one can observe the lack of statistical stability. 
However, this estimate does not show the more irregular behavior during the crisis than during 
the pre-crisis period. 

4. Hill Estimate. (Cf. [5], chapter 3.) If the extreme value index is positive, the Hill estimate 
can be also applied. It is recommended in [3] to test heavy tails using this estimate. For the given 
Xt one can conclude: Xt has a positive extreme value index (within the range 0.3 — 0.8, it slightly 
increased during the crisis). Also, xt has a power tail, i.e. is not normal. In addition, a statistical 
stability is not observed. 



Another drawback of the Samuelson model is that x t are independent. There is an opinion that 
this is not always true. How shall we test whether x t are independent? Let us suggest the following 
test, which we call an M-test. (See also [2].) 

5. M-test. Suppose (X n ,n G Z) is a sequence of independent identically distributed random 
variables with continuous distribution function. Let us call a point n G Z a local maximum point 
if X n > X n _i, X n > X n+ i. Similarly, let us call this point a local minimum point if X n < 
X n +i, X n < X n _\. 

It can be easily shown that local maximum points form an infinite set with infinite upper and 
lower bounds. The same is true for local minimum points. These points alter, i.e. there is a unique 
point of local minimum between any two adjacent points of local maximum, and vice versa. 

Suppose is a local maximum point. Let T\ be the distance from this point to the next local 
minimum point, let r 2 be the distance from this point to the next local maximum point, etc. It is 
straightforward to prove that r = (r ra , n G N) is a strictly stationary sequence and m := Eri = 3/2. 

Conjecture. This sequence is ergodic. 

If this is true then this sequence satisfies the Birkhoff Ergodicity Theorem: 
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almost surely and in the mean sense. The Monte-Carlo modleing showed that this Law of Large 

Numbers is indeed true; in addition, the following limit theorem is true: £ ra := y/n(rh n —m) d] — >Y, 
where EY = 0, VarY ?» 0.73. We have also found some fractiles of Y. It is now obvious how to 
construct an i.i.d. (and continuity of distribution function) test. 

Test results do not entirely conform with the Samuelson model. The dollar has the most irregular 
behavior. It is noteworthy that the results during the crisis do not deviate from m more than the 
results during the pre-crisis period. 

6. Self-similarity Parameter Estimate. The R/S-method is proposed in [3, 4] to estimate 
the self-similarity Hurst parameter for the given process (X n ,n G Z + ) and to test whether this 
process has any self-similarity Let us make its brief outline. Let 
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be an empirical variance for the observations X , . . . , X n _i and let S\. := X/j=o -^-k- The variable 
R{n) = max I Sk S n \ — min I S^ S n 



k=0,...,n-l \ n I k=0,...,n— 1 \ Tl 

is called the range. It is known that the standardized range 
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as n — > oo. Hence one can estimate H by taking logarithms of both parts of this approximate 
equality and applying the method of minimal squares. We have H = 1/2 for X^ := Xk, where Xk 
are from the Samuelson model. 

The following results were obtained: the calculated estimates H contradict the supposition that 
H is equal to 1/2. The estimate is substantially higher during the crisis than for the pre-crisis 
period. (The only exception is the Ukranian griven: H p» 0.927 before the crisis, H p» 0.440 during 
the crisis.) 

All these results clearly show that currency rates dynamics cannot be satisfactorily described 
by the Samuelson model. 

7. Normal Inverse Gaussian Distribution. Using non-Gaussian (e.g., generalized hyperbolic, 
GH) distributions for Xt is a modern approach. Istigecheva (2006, p. 11) claims that these distributions 



can be fruitful in financial time series modeling. They have attractive properties: they are not 
necessarily symmetric and their tails are heavier than normal tail. 

Normal Inverse Gaussian (NIG) distribution is a particular case of GH-distribution. Its density 

is 
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f(x; a, (3, 5, //) := \ . == exp ( 5^a 2 - (3 2 + (3(x - //) ) . 

7T y/8 2 + (X- /J,) 2 V / 

It depends on four parameters: a, (3, 8, \x. Here K\ is a modified Bessel function: 

Kl{x) '-=\L exp ("K y + y)K 

To estimate these parameters, the maximal likelihood estimate is used in the article [1]. But we 
shall use the method of moments. Suppose E, V, S, K is an expectation, variance, asymmetry and 
excess coefficients, then (see [1]) 

5 2 (1 + t 2 ) „ 3r „ 3 / - 2 



E = h + 8t, V= ~ y ~ . ' ' \ S - - = K = - I 1 + 4 ' 



c ' ' Vt^+^Y c V i+T 
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Here we denote /3 = (t/5, a = £vl + t 2 /5. One only needs to calculate estimates for E, V, S, K 
and to solve this system of equations. 

The results: this model is not suitable for currency rates dynamics, since ( has a strong 
statistical instability. 
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